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1. Introduction 

I would like to address two topics in this talk, concerned with Yang-Mills theory and con- 
finement in 2+ 1 dimensions. First, there are now a number of different proposals concerning the 
form of the Yang-Mills vacuum wavefunctional in D = 2 + 1 dimensions, and I will review here 
some recent numerical tests of these proposals. This is work done in collaboration with Hrayr 
Matevosyan, Stefan Olejnfk, Markus Quandt, Hugo Reinhardt, and Adam Szczepaniak, and is pre- 
sented in much more detail in ref. [jl|]. My second topic concerns the Gribov horizon in Coulomb 
gauge. Specifically, I would like to answer this question: Suppose we generate transverse gauge 
fields arbitrarily close to the first Gribov horizon, drawn from a simple momentum-space probabil- 
ity distribution dictated by the desired (e.g. Gribov-Zwanziger) form of the gluon propagator. Are 
such configurations, by virtue of their proximity to the horizon, sufficient to produce a linear (or 
near-linear) rise in the color Coulomb potential? 



2. The Proposed Vacuum Wavefunctionals 

Since most of the interesting properties of non-abelian gauge theories, such as confinement 
and chiral symmetry breaking, are actually properties of the vacuum, we might learn more about 
those phenomena if we knew the explicit form of the ground state wavefunctional [A] ■ A very 
old proposal [0] is that, at large scales, the Yang-Mills vacuum wavefunctional looks like 



% ff [A]=^Vexp 



X -\i Jd 3 xTr[F?j} 



(2.1) 



This vacuum state has the property of "dimensional reduction," in the sense that computation of 
a large spacelike loop in 3 + 1 dimensions reduces to the calculation of a large Wilson loop in 3 

(3) (2) 

Euclidean dimensions. Suppose *J/ () is the ground state of the 3+1 dimensional theory, and *J/q 
is the ground state of the 2+ 1 dimensional theory. If these ground states both have the dimensional 
reduction form, and W(C) is a large planar Wilson loop, then the area law falloff in D = 3 + 1 
dimensions follows from confinement in two Euclidean dimensions in two steps: 

W(C) = (Tr[f/(C)]) D=4 = (^ 3) |Tr[f/(C)]|^ 3) ) 
~ (Tr[U(C)]) D =' = (¥^\Tr[U(C)]\¥^) 

~ (Tr[t/(C)]) D=2 (2.2) 



In D = 2 dimensions the Wilson loop can of course be calculated analytically, and we know there 
is an area-law falloff, with Casimir scaling of the string tensions. 

On the other hand, dimensional reduction cannot be exactly right. For one thing, there is no 
color screening in 2D Yang-Mills, so we would get the wrong N-ality properties in 2+1 dimensions. 
For another, the dimensional reduction state has the wrong high-momentum behavior. Thus there 
must be corrections, and in fact it can be shown [^j], in the context of strong coupling lattice gauge 
theory, that small corrections to the dimensional reduction vacuum wavefunctional are responsible 
for color screening. In this talk I will be concerned with the vacuum state in 2+1 dimensions, where 
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the following well-motivated proposals for the ground state have been advanced by various groups: 



W GO [A] = exp 

*¥kkn[A} = exp 

^hybridiM = exp 

*Pcg[A] = exp 
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(2.3) 
(2.4) 
(2.5) 
(2.6) 



The temporal gauge wavefunctional ^go was suggested by Olejnfk and myself [Qj. In this equation 
— D 2 is minus the covariant Laplacian, Ao is its lowest eigenvalue, and m is a parameter with 
dimensions of mass. Fixing m to obtain the known string tension, we then find that the mass gap 
and Coulomb gauge ghost propagator agree with standard Monte Carlo results. The wavefunctional 
^kkn was obtained by Karabali, Kim, and Nair (KKN) [Q] as an bilinear approximation to the 
ground state within their "new variables" approach. When this approximation is transformed back 
to the usual variables in temporal gauge, the result is the expression shown above. In their case 
m = g 2 N/2n. If one simply drops the term —V 2 , then the string tension agrees remarkably well 
(in the large N limit) with the Monte Carlo result. A questionable feature of this procedure is 
that ^kkn [A] is not gauge-invariant as it stands, and is therefore not a physical state. A reliable 
computation of the string tension via a dimensional reduction approximation really requires, and 
depends on the choice of, a gauge-invariant completion of this state, which goes beyond the bilinear 
result offered by KKN. 1 A gauge-invariant ground state combining features of ^go an d ^kkn 
is the hybrid state ¥ 'hybrid* suggested in [jl|], which agrees with ^kkn when evaluated on abelian 
([Ai,A2] = 0) configurations. Both ^co an d ^hybrid have the dimensional reduction property, when 
evaluated on long-wavelength configurations. The Coulomb-gauge wavefunctional ^cg has been 
advocated and developed by Szczepaniak and co-workers and by Reinhardt and co-workers, cf. 
eg. []|] and [gj. In their proposal co(k) is determined by a set of coupled integral equations. This 
approach leads to an enhancement of the Coulomb ghost propagator, and a confining Coulomb 
potential. On the other hand, an area law for spacelike Wilson loops is not obtained. 

All of these proposals have the same free-field behavior for field configurations with only 
short-wavelength, quasi-abelian ([Ai,A2] ~ 0) components. They differ at the long-wavelength 
end. Therefore we would like to compare the predictions of these proposed wavefunctionals with 
those of the true vacuum wavefunctional [A] evaluated on configurations which are either very 
long wavelength abelian, or essentially non-abelian. 2 



1 For example, a naive replacement of —V 2 by — D 2 in ^kkn would lead to an infinite string tension in the continuum 
limit, cf. 

2 It should be noted that the true vacuum wavefunctional in Coulomb gauge is the same as the true ground state in 
temporal gauge, when both are evaluated on field configurations fixed to Coulomb gauge. 
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3. The Measurement Method 

Consider a modified lattice Monte Carlo procedure in D Euclidean dimensions, in which the 
lattice fields on a particular time slice, t = say, are restricted to be one of a finite set of D — 1 
dimensional lattice configurations = {U^ (x),m = 1,2, ...,M}. Link variables at time t ^ are 
updated normally, but on the time slice I = 0a member of the set is selected at random, and 
accepted or rejected according to the Metropolis algorithm. Let N„ be the number of times that the 
«-th member of the set is accepted. Then it is not hard to show that [|T]] 

where [U ] is the ground state of the lattice transfer matrix in temporal gauge. Thus we are able, 
in principle, to compute the relative amplitudes of the true Yang-Mills vacuum wavefunctional, in 
any given set of configurations. 3 This means that, given any proposal 

x ¥ P ro P o S ai[U]=^e- 1 i R W (3.2) 
for the Yang-Mills vacuum, all we have to do is to plot 



Tog 



vs. RpW] , (3.3) 



where N to t is the total number of updates at t = 0, and serves as a convenient normalization. If the 
proposal is correct, then the data should fall on a straight line with a slope = 1 . 

4. Results 

We will consider three types of trial configurations. The first is abelian plane waves 



U[ m \n x: n 2 ) = ^l-(aW(n 2 )yt 2 + iaW(n 2 )o 3 
U { 2 m \n u n 2 ) = 1 2 

fl (»)( n2 ) = I^TF^cos(^) , (4.1) 

where m = 1 , 2, . . . , m max with L the lattice extension, a , y some constants , and & 2 = 2 ^ 1 — cos ( ^ ^ 
is the squared lattice momentum. To convert to physical units, we will arbitrarily set the string 
tension to a = (440 MeV) 2 , and then the lattice spacing is a = >V G\ at / G as usual. For the GO, 
KKN, and CG wavefunctionals one has 

-logf =2(a + yn)(Q(]?) + ro, (4.2) 



tot 



with the predictions 
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This method was first suggested in ref. [p 
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Figure 1: Cumulative data for (Omc vs - P m physical units, compared to theoretical predictions. 

For CG, a numerical approach is required. In that case the integral equations contain some constant 
parameters c\,C2, and we have w(0) = for c\ = 0, with ft)(0) > otherwise. 



From a best linear fit to the lhs of (42) at a given value of k 2 , we can extract co(k) from the 
data, and compare these values, obtained at various k 2 , to the theoretical predictions. For GO and 
KKN there is a dimensionless parameter g 2 /m which can be chosen from a best fit to the data. Our 
results are shown in Fig. p] for lattices of extensions L = 16,24,32,40,48, and Euclidean lattice 
couplings Pe = 6,9, 12 for the numerical simulations. The GO and KKN (= hybrid for abelian) 
proposals both work well, and are indistinguishable in this range of momenta. The CG proposal 
also works well, for the choice (shown here) of c\ = 0. 

Next we consider non-abelian constant configurations, beginning with configurations in which 
links U\ and U2 have a fixed amplitude, but whose commutator [1/1,1/2] depends on an angle 6, 
which varies in the set: 



u[ m \ni,n 2 ) = \f \ — a 2 tz + iaa\ 

[/^ (711,712) = \/l - a 2 l2 + za(cos(0 m )ai +sin(0 m )a 2 ) , (4.4) 

and we take evenly spaced 8 m = y{m — 1). For the GO and hybrid wavefunctionals, the exponent 
is proportional to the field-strength squared, i.e. R[U] °< (A\ x A2) 2 , for small a, and therefore 

RGO,hybrid[U {n) ]°< sin 2 (9 n ) (4.5) 

for these field configurations. For the Coulomb gauge wavefunctional, however, the exponent is 
proportional to the gauge-field squared, i.e. R[U] °< A 2 +A 2 , and hence, since the amplitudes of A\ 
and A2 are fixed in the set (fO|), 

/? CG [t/W]occo(0), (4.6) 

independent of the angle 8 n . If ft)(0) = 0, which seems optimal for agreement with the plane wave 
data, then Rqg would also be independent of the amplitude of the gauge fields. These predic- 
tions are easy to check, since —\og{N n /N tot ) should be equal to R[U] plus a constant. Plotting 
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—log(N n /N tot ) vs. sin 2 (0„), we find the result shown in Fig. ||. This result is consistent with the 
GO and hybrid proposals, but is clearly not compatible with CG. 



o 




sin 2 (8) 



Figure 2: Dependence of —log(N n /NT) on the "non-abelianicity" of the non-abelian constant configura- 
tions, determined by sin(0„). 



Finally, we consider sets of non-abelian constant configurations with maximal non-abelianicity, 
i.e. 6 = 7r/2, but varying amplitude: 



U{ '{ni,n 2 ) 



l_( a M)2 l2 + /a M ai 



( m ) I \ 
2 («i,"2j 



l_( a M)2 l2 + /a M a2 

n 1/4 

oc + ym 
20L 2 



(4.7) 



We then plot log N m /N tot vs. R[U^ m '], to see whether the data points fall on a straight line and, if 
so, whether the slope equals one. 

An example of the —log[N n /N tot ] vs. RcolU^] data at f$E = 6 is shown in Fig. |3(a)| , for the 
choice a = 2, y = 0.15. Although the data is nicely fit by a straight line which has a slope close to 
unity, this fact must be interpreted with caution because, since the number A^,, falls off exponentially 
with Rgo[U("'], the range of R must necessarily be kept small; typically AR m 4 — 5. This could 
mean that the tendency of the data to lie on a straight line is misleading, and perhaps we are simply 
looking at the tangent of a curve. It is therefore necessary to extract the slope of the straight line 
over small intervals centered around points over a wide range of R. The slope of the data vs. R[U] 
is displayed in Fig. 3(b)| for both the GO and hybrid proposals, and we see that in both cases the 
slope is close to unity. 



5. Numerical Simulation of the GO and hybrid wavefunctionals 

It is also possible to carry out numerical simulations with the probability distribution *}/q [U ] 
for the GO and hybrid proposals, via a method developed in ref. |]3p. Using this method, one can 
compute the mass gap, extracted from the gauge-invariant (F 1 a 2 2 (x)F 1 2 2 (v)) correlator, as well as 
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non-abelian constant, (3=6, L=32 
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Figure 3: (a) Plot of —log(N n /NT) vs. Rqq for non-abelian constant configurations, maximal non- 
abelianicity, at [5e = 6, L = 32, a = 2, y = 0.15. In this case the straight line fit has a slope = 0.98. 
(b) /3 j e= 1 2 calculation, for both types of wavefunctionals. 



gauge-dependent quantities such as the Coulomb gauge ghost propagator 



G(R) 



1 



VD 



(5.1) 



\x-y\=R 



and the color Coulomb potential 

V C (R) 



1 : <-n 1 



V D V V D 



(5.2) 



\x-y\=R 



In the case of the mass gap, the values found from *Pgo appear to converge nicely, in the con- 
tinuum limit, to the values found by standard Monte Carlo techniques in ref. [S], as seen in Fig. 



4(a)[ An even more striking agreement is found between the Coulomb gauge ghost propagator, ex- 
tracted from the probability distributions of *F^ and ^ iyhrid , and the same quantity computed by 
the usual Monte Carlo simulation techniques in three Euclidean dimensions. These are displayed 



in Fig. |4(b)| . Agreement for the color Coulomb potential is not as good, but we have traced the 
disagreement to the fact that this quantity is extremely sensitive to infrequent "exceptional" config- 
urations, in which the value of Ao is unusually low. These rare cases can be eliminated by imposing 
a lower bound cutoff on the value of Vc(0) obtained from a single configuration. As the lower 
bound is raised, agreement between the GO/hybrid and Monte Carlo values is restored, as seen in 
Fig. ||. The color Coulomb potential involves two powers of the inverse Faddeev-Popov operator, 
while the ghost propagator is a single inverse power. We interpret our results to mean that ^ G0 and 
^hybrid are dose to the true ground state ^ 2 true for the bulk of the probability distribution, but may 
deviate somewhat in the tail of the distribution. 



6. Is the Horizon enough? 

I will now turn to the second topic mentioned in the Introduction, relating to the Gribov hori- 
zon. We have seen that the GO and hybrid wavefunctionals provide an area-law falloff for spacelike 
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glueball mass: W G0 vs D=3 Monte Carlo 
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Figure 4: (a) Mass gap extracted from simulations of the GO wavefunctional, compared to the results ob- 
tained by Meyer and Teper [^J. (b)The ghost propagator derived from standard Monte Carlo (MC) simulation 
at Pe = 9, and the same quantity calculated by simulation of the GO and hybrid wavefunctionals. 
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Figure 5: Data for the Coulomb potential at j3g = 9 and L = 32, derived from MC, GO and hybrid simula- 
tions, with a cut on the data, discarding configurations for which |V(0)| is greater than 5, 10, 50, and 300, 
respectively. 



Wilson loops, as well as a confining Coulomb potential, and pass the tests described in previous 
sections. But do we really need such sophisticated wavefunctionals to obtain a confining Coulomb 
potential? Is it not possible that a simple gaussian distribution for momentum-space field com- 
ponents is sufficient for Coulomb confinement, provided the gauge fields are close to the Gribov 
horizon? 

We recall that the Gribov horizon is a cutoff in the functional integral over transverse A-fields, 
which is simply ignored in ordinary perturbation theory. The domain of the functional integral 
should be restricted to configurations for which the Faddeev-Popov operator M = — V • D is positive 
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definite, i.e. the lowest non-trivial eigenvalue Ao of M is positive. This domain is known as the "Gri- 
bov region," and it is bounded by first "Gribov horizon," where Ao = 0. In the Gribov-Zwanziger 
picture of confinement, most of the volume of the Gribov region is concentrated very near the hori- 
zon, and therefore Ao is typically close to zero. Because the Coulomb potential (M _1 (— V 2 )M _1 ) 
contains two inverse powers of the Faddeev-Popov operator, this proximity of typical gauge field 
configurations to the horizon is conjectured to enhance the Coulomb potential from 1/r to linear 
behavior. It was also argued by Gribov that the restriction to the Gribov region would suppress the 
transverse gluon propagator in the infrared, such that in Coulomb gauge, at equal times, 

§. . _ kik i 

Dff(k) = — , 11 != 8 ah , (6.1) 
11 K ' 2^k 2 + m*/k 2 

where m is a constant with dimensions of mass, and k = \k\ refers to the spatial components of 
momentum. 

This leads to an interesting question: Suppose we generate transverse gauge fields which are 
(arbitrarily) close to the horizon, and which result in whatever transverse propagator (e.g. Gri- 
bov's) is desired. Would such gauge fields result in a confining Coulomb potential? To answer 
this question, let us generate some large set of random numbers from a normal distribution, and 
use these random numbers to construct transverse momentum-space gauge fields A,(vt) with the 
required properties. Then we calculate the Coulomb potential, and investigate how that potential 
depends on the proximity to the Gribov horizon, and the form of gluon propagator chosen. 4 

In order to calculate the Coulomb potential, we only need the gauge field at a fixed time. In 
what follows I will use continuum notation, but lattice regularization, and lattice momenta, are 
implicit. From transversality in two dimensions we may write 

A a j(k) = Sj(k)A a (k), (6.2) 

where Ej(k) is the polarization vector, and we want to select A a (k) stochastically from a probability 
distribution such that, on an L x L lattice 

{Al{k)A){k')) = 8 ah ^-^L 2 8 k ^ , (6.3) 

where 

2 \Jk 2 + m 4 /k 2 Gribov propagator 



(0(k) = < 



(6.4) 



+ m 2 massive propagator 
This is achieved by choosing, for each momentum and color component, 



A°(k) = —L=[ n f{k) + i^{k)] , (6.5) 



4 Related work has been carried out in refs. [g| and Q10| ], In those articles the A-fields are drawn from a probability 
distribution corresponding to a Gaussian wavefunctional in 3+1 dimensions, leading to a particular Coulomb gauge 
gluon propagator (infrared finite in the former case, Gribov form in the latter). Transverse gauge fields drawn from 
these distributions were used to calculate the Coulomb gauge ghost propagator numerically, and the results showed at 
best a modest infrared enhancement of the ghost dressing function. Our numerical approach is similar, except that we 
adjust the proximity to the horizon, and calculate the color Coulomb potential (in 2+1 dimensions) rather than the ghost 
propagator. 
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Gribov propagator, L=28 >. =0.033 



0.5 



-0.5 
-1 

-1.5 
-2 

-2.5 
-3 



20 

-20 
-40 
-60 
-80 
-100 
-120 



6 8 10 12 14 16 18 20 
R 



(a) 

Gribov propagator, L=28 X =0.0033 



6 8 10 12 14 16 18 20 
R 



(c) 



Gribov propagator, L=28 X =0.01 

' ' ' -'V'--"'' 



-10 
-12 
-14 



6 8 10 12 14 16 18 20 
R 



(b) 



Gribov propagator, L=28 X o =0. 0033, 0.01 ,0.033 



20 

-20 
-40 
-60 
-80 
-100 
-120 











»•**** 










X =0.0033 . 




X„=0.01 




^0=0.033 



10 
R 



12 14 16 18 



(d) 



Figure 6: Coulomb potential from the Gribov transverse gluon propagator with Ao = (a) 0.033, (b) 0.01, (c) 
0.0033, (d) combined. 

where tj" 2 are random numbers taken from a normal distribution. We then set A"(k) = £j(k)A a (k), 
Fourier transform back to position space, and calculate the Faddeev-Popov matrix (—V • D) c ^. 

The mass parameter m in co(k) can be adjusted to make the lowest (non-trivial) eigenvalue Ao 
of the Faddeev-Popov operator positive, and as close to zero as desired. Then, using the resulting 
gauge field, we calculate the Coulomb potential ( |5.2[ ) at j8 = 4, averaging over many such configu- 
rations. The question is whether this potential is confining for, e.g., a distribution which produces 
a tranverse gluon propagator as proposed by Gribov, given a sufficiently small value for Ao- 

The result for three different choices of Ao = 0.033,0.01,0.0033 are shown in Fig. |6| There 
is no clear evidence of a linear potential developing as Ao decreases. Instead, the main effect of 
decreasing Ao is a drastic increase in the overall magnitude of the color Coulomb potential. We 
can also check whether using a massive transverse gluon propagator, rather than the Gribov form, 
would make any difference. In fact there is no qualitative difference, as seen in Fig. 0. 



It was argued in ref. QlT| ] that confinement calls for not only a small value of Ao, but also an 
enhanced density of near-zero eigenvalues. This may not be a feature of the test configurations we 
have chosen. It is also evident that the A" (k) components generated by our procedure are com- 
pletely uncorrelated, e.g. (A"{qi)A Aq2)A c k (q-})) = 0, which is obviously unrealistic in the infrared 
regime of a non-abelian gauge theory. 

7. Conclusion 

It is found that both the GO and hybrid proposals for the Yang-Mills vacuum wavefunctional, 
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Figure 7: Same as previous figure with a massive transverse gluon propagator, and Ao = (a) 0.033, (b) 0.01, 
(c) 0.0033, (d) combined. 

which both have the property of dimensional reduction, fit the data for abelian plane waves, and 
non-abelian constant lattices, almost perfectly. The mass gap and the Coulomb gauge ghost prop- 
agator also work out well. The Coulomb gauge proposal *Pcg 1S consistent with the abelian plane 
wave measurements, for the choice c\ = 0. It fails on non-abelian constant configurations. 

We conclude that the data supports the conjecture of dimensional reduction in the infrared; i.e. 
long wavelength vacuum fluctuations in 2+1 dimensions, at fixed time, resemble fluctuations in a 
two-dimensional Euclidean theory. A second conclusion is that proximity to the Gribov horizon 
does not, by itself, seem sufficient to produce a confining Coulomb potential, as we have seen in a 
particular example. Field configurations which do give the confining result must evidently satisfy 
some other conditions, apart from proximity to the horizon. 
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